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Proton-nucleus total reaction cross sections in the optical limit Glauber theory: 
Subtle dependence on the equation of state of nuclear matter 
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We calculate the proton-nucleus total reaction cross sections at different energies of incident 
protons within the optical limit approximation of the Glauber theory. The isospin effect has been 
taken into account. The nucleon distribution is obtained in the framework of macroscopic nuclear 
models in a way depending on the equation of state of uniform nuclear matter near the saturation 
density. We find that at an energy of order 40 MeV, the reaction cross section calculated for neutron- 
rich isotopes significantly increases as the parameter L characterizing the density dependence of the 
symmetry energy increases, while at energies of order 300 and 800 MeV, it is almost independent of 
L. This is a feature of the optical limit Glauber theory in which an exponential dependence of the 
reaction cross section on the neutron skin thickness remains when the total proton-neutron cross 
section is small enough. 
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I. INTRODUCTION 

Reactions of unstable neutron-rich nuclei with a pro- 
ton target are of current interest as such reactions can 
act as a major means to probe the matter densities of ex- 
otic nuclei. If one appropriately selects incident energies, 
protons could be more sensitive to neutron distributions 
than to proton distributions of nuclei. The equation of 
state (EOS) of nuclear matter is also of interest since it is 
essential to understand the saturation property of atomic 
nuclei, the structure of neutron stars, and the mechanism 
of stellar collapse. In a previous work [1] we have the- 
oretically connected the EOS of nuclear matter to the 
proton-nucleus elastic differential cross sections. We lim- 
ited ourselves to high energy of order 800 MeV and used 
the optical limit approximation (OLA) of the Glauber 
theory. We find, at large neutron excess, a strong cor- 
relation between the peak angle in the small momentum 
transfer regime and the density symmetry coefficient L, 
which characterizes the density dependence of the sym- 
metry energy and hence controls the size of neutron-rich 
nuclei [2]. At this high energy, we briefly addressed the 
question of how the total reaction cross sections and the 
EOS parameters connect and did not find a strong cor- 
relation as compared with the case of elastic scattering. 

The purpose of the present paper is to study the con- 
nection between the EOS of nuclear matter and the 
proton-nucleus total reaction cross section in more detail 
and at different energies within the OLA of the Glauber 
theory. We are interested in the total reaction cross sec- 
tion because, for neutron-rich isotopes, it is much easier 
to measure experimentally than the elastic differential 
cross section and uncertainties in the measured reaction 
cross section are expected to be sufficiently small to de- 
duce new information about the EOS parameters. 



In the course of this study, however, we happened to 
find some interesting properties of the calculated total 
reaction cross sections, of which the physical manifes- 
tation remains to be examined. We will thus focus on 
these properties from a purely theoretical point of view, 
instead of deducing the EOS parameters from empirical 
data for the total reaction cross sections. 

In Sec. II, we summarize a macroscopic nuclear model 
that allows us to connect the nucleon density distribu- 
tions and the EOS of nuclear matter. Section III is de- 
voted to description of the reaction cross sections within 
the framework of the Glauber theory. In Sec. IV, we cal- 
culate the reaction cross sections for heavy neutron-rich 
nuclei in a way dependent on the EOS of nuclear mat- 
ter, which in turn are analyzed by using the rectangular 
distributions that allow for neutron skin thickness. 



II. FROM EQUATION OF STATE TO NUCLEI 

In this section, we summarize the way we connect the 
nucleon density distributions with the EOS of nuclear 
matter. Generally, the energy per nucleon of nuclear 
matter can be expanded around the saturation point of 
symmetric nuclear matter as 
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Here wq, uq, and Kq are the saturation energy, the sat- 
uration density and the incompressibility of symmetric 
nuclear matter, n is the nucleon density, and a is the 
neutron excess. The parameters L and Sq characterize 
the density dependent symmetry energy coefficient S{n): 
Sq is the symmetry energy coefficient at n = no, and 
L = Suq (dS / dn) n=no is the density symmetry coefficient. 



As shown in Ref. [3 by describing macroscopic nuclear 
properties in a manner that is dependent on these EOS 
parameters, empirical data for masses and radii of stable 
nuclei can provide a strong constraint on the parameters 
wq, no, and 6*0, while leaving Kq and L uncertain. We 
remark that the isoscalar giant monopolc resonance in 
nuclei (e.g., Ref. [3|) and caloric curves in nuclear colli- 
sions (e.g., Ref. [J]) can constrain Kq only in a way that 
is dependent on models for the effective nucleon-nucleon 
interaction. 

The incompressibility Kq and the density symmetry 
coefficient L control in which direction the saturation 
point moves on the density versus energy plane, as the 
neutron excess increases from zero. This feature can be 
found from the fact that up to second order in a, the 
saturation energy Wg and density n^ are given by 
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are the potential energy densities for symmetric nuclear 
matter and pure neutron matter. This expression for w 
is one of the simplest parametrization that reduces to Eq. 
([1]) in the simultaneous limit of n — >■ no and a — > 0. We 
then set the nucleon distributions n,j(r) (i — n,p) as 
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(3) and in the spirit of the Thomas-Fermi approximation 
minimize the total energy E with respect to i?i, ii, 



The influence of Kq and L on neutron star structure can 
be significant [5] despite the fact that these parameters 
characterize the EOS near normal nuclear density and 
proton fraction, which are fairly low and large, respec- 
tively, as compared with the typical densities and proton 
fractions in the central region of the star. 

We proceed to macroscopic nuclear models used in this 
and previous works [ll, Hi Q • We describe a spherical nu- 
cleus of proton number Z and mass number A within 
the framework of a simplified version of the extended 
Thomas- Fermi theory f7|. We first write the total en- 
ergy of a nucleus as a function of the neutron and proton 
density distributions nn{r) and np(r) in the form 



E = Eb + Eg+Ec + {A- Z)m„c'' + Zmpc"" 



where 



Eb= drn{r)w {nn{r),np{r)) 
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(5) 



is the bulk energy with the energy per nucleon w{nn,np) 
of uniform nuclear matter. 
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is the gradient energy with adjustable constant Fq, 
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is the Coulomb energy, and 7n„ and mp are the neutron 
and proton rest masses. We express w as Q 
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and n™ with the mass number A^ the EOS parameters 
{uq^wq, SqjKqjL) and the gradient coefficient Fq fixed. 
By calculating the charge number, mass excess, and root- 
mean-square charge radius from the minimizing values of 
Ri, ti, and n™ and fitting the results to the empirical 
values for stable nuclei (25 < A < 245) on the smoothed 
beta stability line, we finally obtain no, wq, Sq, and Fq 
for various sets of L and Kq ranging MeV < L < 175 
MeV and 180 MeV < Kq < 360 MeV. We remark that as 
a result of this fitting, the parameters ai-62 become func- 
tions of Kq and L, while we fix the remaining parameter 
63, which controls the EOS of matter for large neutron 
excess and high density, at 1.58632 fm'^ throughout this 
fitting process. 

The macroscopic nuclear models used here can describe 
gross nuclear properties such as masses and root-mcan- 
square radii in a manner that is dependent on the EOS 
parameters, L and Kq. Notably, as in Ref. [2;], these 
models predict that the root-mean-square matter radii 
increase appreciably with L, while being almost inde- 
pendent of Kq, and that the root-mean-square charge 
radii are almost independent of both L and Kq. How- 
ever, there are some limitations in the present macro- 
scopic approach. First, this approach works well in the 
range of a ^ 0.3 and A <; 50, where a macroscopic view 
of the system is relevant. Second, the nuclear surface is 
not satisfactory in the present Thomas-Fermi-type the- 
ory, which tends to underestimate the surface diffuseness 
and does not allow for the tails of the nucleon distribu- 
tions. Lastly, no shell and pairing effects are included. 



III. GLAUBER MODEL 

In the OLA of the Glauber model [g|, the proton- 
nucleus scattering phase-shift function is simply de- 



scribed by the proton density np{r) and the neutron den- 
sity 7T,„(r) as follows: 



a«XOLA(b) 



exp[iXp(b) +'i'Xnib)], 



(12) 



where Xp (Xn) is the phase shift due to protons (neutrons) 
inside the nucleus, 

iXp{b) = - drnp{r)Tpp{b~ s), 

iXnib) == - j drnn{r)Tpn{b- s), (13) 

with the impact parameter b and the projection s of the 
coordinate r on a plane perpendicular to the incident 
proton momentum. 

The profile function, FpAr, for pp and pn scatterings, is 
usually parametrized in the form: 



^pN{b) = 
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where apN is the ratio of the real to the imaginary part 
of the pp (pn) scattering amplitude in the forward direc- 
tion, cr''^ is the pp (pn) total cross section, and PpN is 
the slope parameter of the pp (pn) elastic scattering dif- 
ferential cross section. These parameters are tabulated 
for different incident energies in Ref. [9J. 

The total reaction cross section of the proton-nucleus 
collision is calculated from 
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IV. 



EOS DEPENDENCE OF REACTION CROSS 
SECTIONS 



In our calculation we have used a code based on the 
Monte Carlo integration for evaluations of the phase shift 
function, which can be applied to any arbitrary form of 
nucleon distributions. We have calculated the total re- 
action cross section for p-^^Cu, p-^"Ni, p-^^^Sn, p-^^'^Sn, 
and p-^°*Pb by substituting the distributions (dJ) deter- 
mined as functions of L and Kq into n„ p(r) in Eq. (|13p. 
For each system at each energy we have performed sum- 
mation over 228 random points, which ensures sufficient 
accuracy. 

Figure [T] shows the results for the selected isotopes at 
800 MeV. As we pointed out in our previous paper [l\ 
the dependence of the total reaction cross section on the 
EOS parameters is rather weak. By decreasing the en- 
ergy to 300 MeV, we obtain the same conclusion as at 
800 MeV; the results are shown in Fig. [2] Note that 
at such high energies the reaction cross section does not 
obey an empirical law obtained within a black sphere ap- 
proximation [10|, which suggests that the reaction cross 
section increases with the nuclear size and thus L. This 
is a feature to be discussed in more details below. 
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FIG. 1: (Color online) The total reaction cross sections cal- 
culated as a function of A'o and L for p-^^^'^^^Sn and p-^'^Ni 
at 800 MeV. 



At 40 MeV the cri°' and cri°' are fairly large compared 



^ pn 



PP 



to those at 800 MeV, and cr*°' is much larger than 0-*°'. 
Therefore we expect for neutron rich nuclei that the total 
reaction cross section depends sensitively on the neutron 
density distribution including the skin region. 

Figure [3] shows our calculations of the reaction cross 
section as a function of L and Kq. We find, for neutron 
rich isotopes like ^"Ni {Z/A = 0.35), a strong depen- 
dence on the L value. The value of the total reaction 
cross section for {L,Ko)=il2A, 350) in MeV is 1248 mb, 
while that for {L,Ko)={171.3, 360) in MeV is 1346 mb. 
The difference between them is about 7%. If uncertain- 
ties of experimental data for the reaction cross section are 
less than about 3.5% and our calculations with simplified 
scattering theory and density distributions are satisfac- 
tory, we could determine the L value from comparison 
between the calculated and measured values. 

For possible determination of L, therefore, it is impor- 
tant to examine how well the OLA of the Glauber model 
predicts the reaction cross section. Although this de- 
pends on the choice of density distributions, for simplic- 
ity, we here confine ourselves to rectangular nucleon dis- 
tributions, for which one can calculate the reaction cross 
section analytically in the zero range limit {PpN — ^ 0). 

Let us set the radius and density of the proton (neu- 
tron) distribution for a nucleus of given A and Z to be 
o-p (an) and n^o — iZ/Airap [n„o = 3(A — Z)/47ra^]. If 
Q-B = On = o,, one can calculate the reaction cross section 



15)1 in the zero range limit as 
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(16) 



where ^ = 2(cr*°'npo + o-p°'n„o)a is the effective optical 
depth. Note that in this case, the OLA of the Glauber 
model is strictly applicable when the effective optical 
depth is close to zero and the proton wave length is 
sufhciently short. In the limit of complete absorption 
(^ — >■ oo), however, Eq. (J16p reproduces the correct an- 
swer 7ra^ . One may thus expect that the OLA provides a 
fairly good prediction of the reaction cross section even 
for intermediate values of ^ as long as the incident energy 
is sufficiently high. Nevertheless, Eq. P^ . if the expo- 
nential term in the right side works, is at odds with the 
power-law ct*'JJ dependence of the reaction cross section 
that is empirically suggested [ll| . Even if the exponential 
term is negligible, the energy dependence of the reaction 
cross section as deduced from Eq. ([T6| shows only a weak 
dependence on A, which is again at odds with the empir- 
ical behavior [12l]. 

We now consider the case in which a„ and ap are gen- 
erally different from each other. In this case, as we shall 
see, the reaction cross section (jlSp in the zero range limit 
includes a term depending exponentially on the neutron 
skin thickness a„ — ap. In fact, when one retains correc- 
tions to the black-disk limit {crl^^ — > oo) by transparency 
of the skin region, the reaction cross section (|15p in the 
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FIG. 3: (Color online) The total reaction cross sections cal- 
culated as a function of Kg and L for p-^^^'^^^Sn and p-^^Ni 
at 40 MeV. The experimental values (solid lines: central val- 
ues, dashed lines: upper and lower bounds) are taken from 
Ref. Q. 
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FIG. 2: (Color online) The total reaction cross sections cal- 
culated as a function of Kq and L for p-^^^'^^*Sn and p-*'^Ni 
at 300 MeV. 
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2crp° n„o- 



where C„ 

It is instructive to note that Eq. ([T7|) reduces to 7ra^ 
in the black-disk limit. For large but finite values of C„, 
the remaining term in the right side of Eq. (J17p starts 
to play a role ahead of the corrections responsible for 
transparency of the inner region in which protons are 
present. Interestingly, this term is a decreasing function 
of a„ and in some cases acts to cancel an increase of the 
main term tto^ by an increment of a„. As we shall see, 
this term is responsible for the L dependence of ctr as 
shown in Figs. 1-3. 

Let us now assume that the distribution (ITT|) roughly 
corresponds to a rectangular one discussed above. Then, 
one can set the L and Kq dependence of a„ and ap in 
such a way that o„ increases linearly with L, while any 
other dependence is negligibly weak (see Fig. 7 in Ref. 
[3]). For ^''Ni, in particular, we set a„ and ap as 
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FIG. 4: (Color online) The total reaction cross sections cal- 
culated from Eq. ^ for '^"Ni by using Eqs. gill and ^ 
as a„ and ap. The full, long-dashed, and short-dashed lines 
are the results obtained for o-p^* = 220, 35, 38 mb, which are 
relevant at incident energies of 40, 300, 800 MeV [11]. For 
comparison, we also plot by a dotted line the result obtained 
for (j*°' = 10 mb. 
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OLA calculations. 

In fact, as shown in Fig. 4, we find that Eq. ([T7|) for 
various sets of ap'^ gives a roughly linear dependence 
on L, while this dependence is suppressed for small val- 
ues of CT*°* . This is consistent with the tendency seen in 
Figs. 1-3, although higher order corrections to the black- 
disk limit have to be significant for small values of Cn or, 
equivalently, small cTp'll^. 

Figures 1-4 suggest that the magnitude of the reaction 
cross sections depends on the incident energy more re- 
markably in the OLA calculations than in Eq. (|T7l) . This 
is because in the OLA calculations we used the density 
distributions of form ([TT|) . which has a surface diffuse- 
ness and thus allow reactions to occur in an even farther 
region than the radius of the rectangular distributions. 

In summary, we have discovered that the OLA of the 
Glauber theory exhibits a nonnegligible exponential de- 
pendence of the total proton-nucleus reaction cross sec- 
tions on the neutron skin thickness at incident energies 
where (7*°* is sufficiently small. It is, however, important 
to note that the validity of using the OLA of the Glauber 
theory in evaluating the total reaction cross sections in 
the energy range considered here is far from obvious. We 
also note that the Coulomb force ignored here has to play 
a role in distorting the proton trajectory at the lowest en- 
ergy. For duly describing the energy and size dependence 
of the total reaction cross sections, therefore, alternative 
approaches based on empirical data for the total reac- 
tion cross sections such as those in Refs. [ll|, [l5| might 
be useful. 



in reference to the root-mean-square neutron and proton 
radii calculated in Ref. 2']. Here we have taken into ac- 
count a factor -\/5/3, which is the ratio of the half-density 
and root-mean-square radii for rectangular distributions. 
By using such a rough correspondence, we can under- 
stand the peculiar L dependence of ctr obtained in the 
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